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ABSTRACT 

The two edge-on galaxies NGC 891 and NGC 7814 are representative of two extreme morphologies: the former is disk-dominated 
while the latter is almost entirely bulge-dominated. It has been argued (van der Kruit 1983) that since the two galaxies, which are 
optically so different, have similar rotation curves their total mass distributions cannot be related in any way to the light distributions. 
This would lead to the conclusion that dark matter is the dominating component of the mass. We have derived new rotation curves 
from recent, high-sensitivity H I observations and have found that the shapes of the rotation curves are significantly different for the 
two galaxies. They indicate that in NGC 7814 the mass is more concentrated to the centre as compared to NGC 891. This reflects the 
distribution of light which is more centrally concentrated in NGC 7814 than in NGC 891. Mass and light do seem to be closely related. 
This is confirmed by the analysis of the rotation curves in mass components: solutions close to the maximum light (bulge + disk) do 
provide excellent fits. In NGC 891 bulge and disk can explain the rotation curve without any need for dark matter out to ~ 15kpc. 
In NGC 7814 the bulge dominates in the inner parts; further out the rotation curve is well reproduced by a maximum disk but its 
M/L ratio is excessively high. A substantial dark matter contribution, closely coupled to the luminous component, seems, therefore, 
necessary. 

Key words, galaxies: individual (NGC 891, NGC 7814) - galaxies: kinematics and dynamics - galaxies: structure 



1. Introduction 

The distribution and relative importance of luminous and dark 
matter in galaxies are still a matter of debate. The distribution of 
mass in a spiral galaxy is inferred from its rotation curve. In the 
past three decades, rotation curves have been derived for a num- 
ber of spirals of various masses and morphological types. For the 
interpretation of the rotation curves in terms of mass distribution 
use has been made of multicolour photometry. The conclusion 
has been that in the outer parts of galaxies there is a significant 
discrepancy between the observed curve and the curve predicted 
from the photometric and gas profiles. Such discrepancy is usu- 
ally interpreted as evidence for the presence of dark halos around 
spiral galaxies. For the bright inner parts of the disk, inside R25, 
there is no such con sensus. There is a range of possibilities from 
"maximum disks" dKalnai si 1 1983b Ivan Albada & SancisiU l986: 
lKentL[l98r3) with constant values of the M/L ratio (implying that 
luminous matter dominates) to "minimum disks" with dark mat- 
ter being the dominant component everywhere. 

The debate on the relative distribution of luminous and dark 
matter in spiral galaxies is still open. The comparison between 
the two edge-on galaxies NGC 891 and NGC 7814 provides a 
good illustration. These two spirals are representative of two ex- 
treme morphologies: a disk-dominated NGC 891 and an almost 
entirely bulge-dominated NGC 78 14 (see Fig. [TJ. Van der Kruit 
(1983, 1987, 1995) pointed out that, in spite of this striking dif- 
ference in their light distributions, these two galaxies have, be- 
yond the central 3 -4 kpc, "essentially identical" rotation curves. 
These curves were measured in the H 1 21 -cm line but they were 
not known inside 1 arcmin (3 - 4 kpc) because HI emission 
was not detected there. The conclusion was that "if the distri- 
bution of luminous matter were in any way related to the total 



mass distribution it would not be possible for the disk-dominated 
and the bulge-dominated galaxies to have such similar rotation 
curves". If the two galaxies had indeed rotation curves of such 
similar shape this would be contrary to the rule that in spiral 
galaxies there is a close c orrelation betw een light distribution 
and rotation curve shape (ISancisiL 12004]) . The contrast would 
also exist with non-Newtonian theories o f gravity which woul d 
predict grossly different rotation curves (Ivan derKruitl |l995). 
There is little doubt that the light distributions of NGC 891 and 
NGC 7814 are very different. But, do their rotation curves have 
indeed essentially identical shapes? 



The rotation curves used by 



Van 
and 



der Kruit for 
for NGC 78 14 



NGC 891 dSancisi & Alien! \l9ld 
(Ivan der Kruit & Searle, Il982l) were obtained from H I ob- 
servations with the Westerbork Synthesis Radio Telescope 
(WSRT). The rotation curve of NGC 891 was obtained directly 
from the observations using a velocity "envelope" method, that 
of NGC 7814 by constructing models. Both curves appear to be 
flat to a first approximation, but the H I observations, especially 
for NGC 7814, had poor signal/noise ratios in the central regions 
and the curves had large uncertainties or could not be derived at 
all. 

In recent years, new H I observations were obtained for these 
two galaxies with better sensitivity with the WSRT. From these 
observations we derived new rotation curves. We used Spitzer 
3.6 /urn observations to derive the photometric parameters of the 
bulge and the disk components and perform mass decomposi- 
tions. We were able to reach new conclusions on the relative dis- 
tribution of luminous and dark matter, which differ from those 
mentioned above. 
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Fig. 1. Left: Multi-band (Y, J, K) image of NGC 891 obtained with the WIRCam at the Canada-France-Hawaii Telescope (CFHT). 
Right: Multi-band (g, r, i) image of NGC 78 14 from the Sloan Digital Sky Survey. 

The derivation of the rotation curves was done in two steps. 
A first estimate was obtained by taking the envelope on the high 
rotational velocity side of the position-velocity diagram along 
the major axis. For a descriptio n of the method and a d iscussion 
of the uncertanties involved see Sancis i & Alien! d 19791) . The ro- 
tational velocity was derived by fitting Gaussian functions with 
fixed dispersions to the high rotational velocity sides of the line 
profiles (i.e. the lowest radial velocities on the approaching side 
and the highest on the receding one). It was assumed that the 
gas is in circular motion, that there is gas emission at the line 
of nodes and that the velocity dispersion is constant and equal 
to 8kms~'. Subsequently this rotation curve was used, together 
with the observed radial H I density profile, as input for the mod- 
elling of the data cube. This was done by assuming concentric 
rings. Each ring has its own rotational velocity and gas den- 
sity. The inclination angle is fixed to 90°, determined by mod- 
elling the HI datacube. The position angle was measured us- 
ing the total HI map and the Spitzer data, the values agreed 
within < 0.5°. There is no evidence for a change of more than 
2° in position angle and 4° inclination out to the outer radius 
from which the rotation curve is deriv ed (for NGC 891 see also 
lOosterloo. Fraternali & Sancisill2007l) . The effect on the rotation 
curves of these changes would be less than 1 kms -1 . 

The gas density as a function of R was obtained by depro- 
jecting the total HI maps. The model position-velocity diagram 
was compared (by eye) with the observed one. When necessary, 
the initial input values of the rotational velocity were changed 
and the procedure was iterated until a satisfactory matching of 
model and observations was reached. This was done indepen- 
dently for the two sides of the galaxies. The final adopted ro- 
tation curves are shown plotted on the position-velocity maps 
(Figs.[2]&[3]). The averaged curves are shown in Fig. For each 
point, the error is the larger between the formal error of the fit 
and an uncertainty due to asymmetries between the approach- 
ing and the receding sides. The latter was calculated assuming 



Table 1. Optical and H I parameters for the two galaxies 



Parameter 


NGC 891 


NGC 7814 


Ref. 


Morphological type 


Sb/SBb 


Sab 


1,2, 


Centre (a 12000) 


2 h 22 m 33.41 s 


h 03 m 14.89 s 




(<5 12000) 


42°20'56.9" 


16°08'43.5" 




Distance (Mpc) 


9.5 


14.6 


1" 


Lb (Lg,b) 


2.5 x 10 10 


1.3 x 10 10 


3" 


Lk (Lo,b) 


7.8 x 10 10 


6.4 x 10 10 


4 


R25 in B-band (kpc) 


18.6 


9.6 


3,5 


R20 in K-band (kpc) 


10.4 


7.2 


4 


Systemic velocity (km s~') 


528 ±2 


1043 ± 4 


6,7 


Total HI mass (M ) 


4.1 x W 9 


1.1 x 10 9 


6,7 


H I inclination (°) 


> 89 


~ 90 


6,7 


1 arcmin = (kpc) 


2.76 


4.25 


7 



(1) 

(3) 



van der Kruit & Searld dl98lh . ( 2)lGarcia-Burillo & Guelii] i 19951) . 
de Vaucouleurs et aLU^giTTwlJarrett et all (120031). (5) 



lAbazaiian et al.ld2004l) . (6l lOosterloo. Frate rnali & S ancisil d2007h . (7) 
this work. 

" More recent determinations of the distances of NGC 891 and 
NGC 7814 using the tip of the red giant branch are fully consistent 
with these values, i.e. 9. 1 ± 0.4 Mpc and 14.8^ '1 Mpc respectively 
dRadburn-Smith et aUl201lh . 
Corrected for internal and Milky Way extinction. 



2. Derivation of the rotation curves 

For NGC 89 1 we used t he H I observations by 
lOosterloo. Fraternali & Sancisil (120071) and for NGC 7814 

(12008). The latter consist of 



those obtained bv iKam nhuis 



4x12 hours integration with the WSRT (September 2004). The 
reduction was done with the MIRIAD package and standard 
calibration was applied. The result is a data cube consisting of 
160 channels with a spatial resolution (FWHM) of 54.5" x 12.9" 
and with a velocity spacing of 4.12 km s . 
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Fig. 2. Upper: Position-velocity diagram along the major axis 
of the edge-on galaxy NGC 891. The squares show the rotation 
curves derived separately for the approaching and the receding 
side. The contour levels are: 0.2, 0.45, 1.0, 2.0, 4.0, 10, 20, 40 
mJybearrT 1 . Lower: Position-velocity diagram along the major 
axis for the model galaxy. The model has been constructed for 
the symmetric part of the disk. 1 arcmin= 2.76 kpc. 



that the d ifference betwe en the two sides corresponds to a 2<x 
deviation dSwaterslll999l) . 

The HI disk of NGC 891 is k nown to be lopsided 
(Bald win. Lvnden-Bell & SancisiL 1980). On its southern reced- 
ing side, it has an extended tail (Fig. [2j at a lower radial ve- 
locity with respect to the central part of the disk. However, 
this decrease in radial velocity does not necessa rily mean a de- 
cline in rotational velocity. As pointed out by ISancisi & Alien! 
(fl979l) it could be due to the location of the gas away from 
the line of nodes. Also, it is possible that non-circular mo- 
tions dominate in these outer parts where the galaxy becomes 
so asymmetric. We derived, therefore, the rotation curve only 
for the symmetric part, inside ~ 6 arcmin (~ 17 kpc) from 
the centre (Table |2). The sharp rise and peak near the cen- 
tre (R < 1 kpc) indicate the presence of a fast rotating inner 
HI disk or ring. Although this seems the straightforward ex- 
planation the possibility of a bar and associated non-circula r 
motions cannot be ruled out dGarcia-Burillo & Guelini 119951) . 
However, in the present data there is no indication of non- 
circular motions. Just outside this central disk or ring the ro- 
tation curve is not well defined due to some beam-smearing 
and blending with the emission from the inner disk. The sec- 



Table 2. Rotation curve of NGC 891 



Radius 


Radius 






(arcmin) 


(kpc) 


( km s 


- 1 \ 

) 


0.32 


0.88 


234.5 ± 


12.0 


0.80 


2.22 


191.8 ± 


12.0 


1.13 


3.11 


211.7 ± 


6.3 


1.45 


4.00 


223.3 ± 


4.4 


1.77 


4.89 


222.5 ± 


4.0 


2.09 


5.78 


223.6 ± 


3.8 


2.41 


6.67 


224.2 ± 


3.2 


2.73 


7.56 


226.1 ± 


3.4 


3.06 


8.44 


226.4 ± 


3.2 


3.38 


9.33 


226.6 ± 


2.7 


3.70 


10.22 


226.6 ± 


2.7 


4.02 


11.11 


223.9 ± 


2.4 


4.34 


12.00 


220.1 ± 


2.4 


4.66 


12.89 


217.6 ± 


2.3 


4.99 


13.78 


217.2 ± 


2.8 


5.31 


14.67 


215.6 ± 


4.4 


5.63 


15.56 


210.1 ± 


4.5 


5.95 


16.44 


207.6 ± 


4.4 



ond point shown here is, therefore, rather uncertain. We used 
different density distributions for the inner ring and estimate an 
error of about 12 km s . In these central parts, the rotation curve 
is slightly asymmetrical: on the receding side it is somewhat 
steeper. Further out the approaching and receding rotational ve- 
locities are the same within < 10 km s . The rotation cu r ve de- 
rived here differs from that obtained by ISancisi & Allenl (fl979) 
especiall y in the inner parts where the higher sen sitivity observa- 
tions bv lOosterloo. Fraternali & Sancisil (120071) used here show 
the inner steep rise and the fast-rotating component. The radial 
extent is about the same. 

The rotation curve of NGC 78 14 (Table [3} has a very steep 
rise near the centre followed by a slow decline which becomes 
a little more pronounced beyond 2 arcmin (~ 8 kpc). The curve 
in the inner 3 kpc from the centre is not accurately determined 
because of the poor signal/noise ratio. We estimate an error of 
about 15kms . The value of the rotational velocity near the 
centre (first two points) could be higher than given here. The 
third point of the rotation curve is well determined (see Fig.0. 
There are no large asymmetries; only at large distances from the 
centre the two sides of the galaxy are not completely symmetri- 
cal. Fig.|3]shows that there are differences in the H I density dis- 
tribution and also in the kinematics. On the approaching side the 
rotational velocity begins to drop off around 2 arcmin from the 
centre whereas on the receding side the decrease seems to start 
a little further out. This rotation curve is different from the flat 
curve derived by Van derKmit andSearle (1982). Because of the 
better signal/noise ratio of the new observations it was possible 
to derive it also in the inner region (within 1 arcmin ~ 4 kpc from 
the centre) where it has a steep rise to a maximum of 250 km s _1 
at about 0.6 kpc from the centre followed by a slow decline to 
about 214kms~ I in the outer parts. Its radial extent is about the 
same. 

3. Comparison of the rotation curves with the 
distribution of light 

The rotation curves for the two galaxies are compared in Fig. [4] 
They have similar amplitudes but significantly different shapes. 
The rotation curve of NGC 7814 indicates a more pronounced 
mass concentration to the centre than in NGC 891. The corre- 
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4. The rotation curves of NGC 89 1 and NGC 78 14. The open symbol for NGC 89 1 outlines some potential contribution from 
■circular motions in the centre. 



Table 3. Rotation curve of NGC 78 14 



Radius 


Radius 


v c 




(arcmin) 


(kpc) 


(kms 


-') 


0.15 


0.64 


250.0 ± 


15.0 


0.41 


1.74 


240.0 ± 


10.0 


0.67 


2.83 


230.6 ± 


7.3 


0.93 


3.96 


230.0 ± 


6.5 


1.20 


5.10 


228.1 ± 


5.4 


1.47 


6.23 


227.9 ± 


5.2 


1.73 


7.36 


226.9 ± 


5.9 


2.00 


8.49 


226.1 ± 


3.0 


2.27 


9.63 


222.8 ± 


3.9 


2.53 


10.76 


218.6 ± 


2.9 


2.80 


11.89 


216.2 ± 


2.6 


3.07 


13.02 


214.2 ± 


2.6 


3.33 


14.16 


213.9 ± 


3.1 


3.60 


15.29 


213.9 ± 


3.6 


3.87 


16.42 


213.5 ± 


4.5 


4.13 


17.55 


213.8 ± 


5.1 


4.40 


18.69 


213.6 ± 


4.8 


4.67 


19.82 


214.3 ± 


4.9 



spondence with the central concentration of light in NGC 7814 
seems evident. Given the potential uncertainties in the interpreta- 
tion of the first point for NGC 89 1 , we indicate it as an open sym- 
bol. In order to further investigate the distribution of the mass as 
compared to that of the light in the two galaxies we made the 
standard decomposition of the rotation curves in mass compo- 
nents: bulge, stellar and gaseous disks, and dark matter halo. 

3.1. Photometric data 

In order to avoid problems with dust extinction we used data 
in the 3.6;um band obtained with the Spitzer Space Telescope 
dWerner et all l2004l) . Both NGC 891 and NGC 7814 had been 
observed with the Spitzer Telescope and the mosaics were al- 
ready available in the archive. The observations of NGC 891 
had an exposure time of 96.8 s, those of NGC 7814 26.8 s. We 
converted from the Spitzer units of MJysr -1 to magarcsec -2 
in 3.6 //m band assuming an a bsolute magnitude for the Sun of 
Mxetim = 3.24 (lQhet all l2008h . We applied the surface bright- 



ness correction simply by multiplying our fluxes by 0.91 (Spitzer 
Handbook). We estimated the values of the background in the 
two images to be 21.06 and 20.70 mag arcsec -2 respectively for 
NGC 891 and NGC 7814. The 3.6 fim images for the two galax- 
ies are shown in Fig. 5 and in Fig. 6. 

We perform ed a bulge-disk decomp osition of the stellar light 
using GALFIT dPeng et al.L 120021 |2Q10|> . For both galaxies we 
extracted the psf directly from the image taking a relatively iso- 
lated star in the field close to the galaxies. The images were 
masked in order to exclude stars, background galaxies and spu- 
rious features of the CCD. We used a Sersic profile for the bulge 
component and a perfectly edge-on disk for the stellar disk. 
NGC 891 could not be fitted with a single disk but it required 
both a thin and a thick disk. The p resence of two disks had been 
found also before in other bands (Ivan der Kruit & Searlel 1.1981; 
IShaw & Gilmorel Il989t IXilouris et al.L 119981) . We note that the 
functional form of the density distribution in the vertical direc- 
tion has a strong impact on the properties of the disk(s) required 
by the fit. At the moment GALFIT allows only for a sech 1 ver- 
tical density profile, but it is possible that with an exponential 
profile the two components could have very different parameters 
or the second disk might even not be required. 

In fitting the bulge and the disk(s) we first fixed the centres 
and the position angles of the various components to the same 
value. For NGC 78 14, we fixed the Sersic index to 4, i.e. we 
assumed a de Vaucouleurs profile; leaving the Sersic index free 
gives values very close to 4. In NGC 891 there is a degeneracy 
between the scale-length of the bulge and the scale-height of the 
thick disk and, to obtain a satisfactory result, we preferred to 
fix the latter to a value of h z = 0.8 kpc. Once these values are 
fixed, the simultaneous fits of all the other parameters converge 
for both galaxies. 

Figs.[5]and[6]show our best photometric models for NGC 89 1 
and NGC 7814 respectively. In both figures the top panels show 
the light distribution (at 3.6 /mi) in color shade and thin (blue) 
contours from the Spitzer mosaic. The thick (black) contours 
show the GALFIT solutions obtained using two disks + a Sersic 
(n=3) bulge for NGC 891 and one disk + a r~ 1/4 bulge for 
NGC 78 14 (see Table g}. The bottom panels in Figs. and 
show ID cuts of the data and the models along the vertical di- 
rections (left) and along the plane of the galaxies (right). In the 
plots perpendicular to the plane we show the data from both sides 
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Fig. 5. 3 .6 fim photometry for NGC 891. Top panel: distribution of light in the sky in thin (blue) contours overlaid with our best 
model with bulge and thin and thick disks, see Table [4] Bottom left: Light distribution along the vertical direction above and below 
the plane (circles and triangles) compared with our best model (thick blue curve). Bottom right: Light distribution along the plane 
of the galaxy (circles and triangles) compared with our best model (thick blue curve). 



(circles and triangles), our best model (thick blue curve) and 
the contribution of the disk only (thin black curve). Note that 
this latter is completely negligible for NGC 7814. In the plots 
along the plane for NGC 891 one can appreciate that the disk in 
not exponential. This shape has been modelled using the trunca- 
tion function imple mented in GAL FIT with rb rea k — 120" and 
Ar so f t = 280" (see lPeng et aUl20Tol) . The truncation is the same 
for both the thin and the thick disk. The disk of NGC 78 14 does 
not require any truncation. 

Table|4]gives the photometric parameters of the fits shown in 
Figs. [5] and [6] Since our rotation curve fitting routine requires 
a spherical bulge we used a spherical-equivalent effective ra- 
dius r e , S ph = r e yfq (geometrical mean) where q is the axis ra- 
tio. The 3.6 /im-band disk and bulge luminosities in Tableware 
calculated out to the last measured points of the rotation curves, 
R m . dX = 17kpc and 20kpc for NGC 891 and NGC7814 respec- 
tively. The values derived for the NGC 891 disks do not take 
into account the truncations and are therefore upper limits. As 
expected, the dominant luminosity component in NGC 891 is 



Table 4. Spitzer 3.6/im-band photometric parameters for the 
two galaxies 



Parameter 


NGC 891 


NGC 78 14 




Thin disk Thick disk 




h R (kpc) 


4.18 5.13 


4.26 


£ E (L 0i 3.6 (1 mPC ") 


443.5 368.9 


78.0 


h z (kpc) 


0.25 0.80° 


0.44 


Alisk ( L 0- 3.6fim) 


4.4 x lO 10 * 5.1X10 106 


8.4 x 10 9 


Sersic index 


2.99 


4.0" 


r e (kpc) 


1.80 


2.16 


h (L Q ,3.6jim P c 2 ) 


525.8 


1.12x 10 3 


q 


0.68 


0.61 




2.2 x 10 10 


7.0 x 10'° 



" fixed in the fit; * calculated out to R = 17 kpc without considering the 
truncation. 
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Fig. 6. 3.6 pm photometry for NGC 7814. Top panel: distribution of light in the sky in thin (blue) contours overlaid with our best 
model with bulge and disk, see Table [4] Bottom left: Light distribution along the vertical direction above and below the plane 
(circles and triangles) compared with our best model (thick blue curve). Bottom right: Light distribution along the plane of the 
galaxy (circles and triangles) compared with our best model (thick blue curve). 



that of the disks, about four times brighter than the bulge. On 
the contrary, in NGC 78 14 the bulge is totally dominating with 
a bulge-to-disk ratio in the 3.6A<m-band of about 9. In K -band 
this ratio is about 14 dWainscoat. Hvland & Freeman! 1 19901) . 

3.2. Maximum light (bulge + disk) 

Fig. [7] (upper panels) shows the 3.6 yum-band photometric pro- 
files adopted for bulge and disk. For NGC 891 the thin and 
thick disks are combined in a single component with inter- 
mediate scale-height. The truncation in GALFIT is applied to 
the projected edge-on profile and it requires deprojection to 
a face-on view. If the parameter Ar so f t is small compared to 
'"break this can be achieved with good accuracy by simply ap- 
plying the truncation function to the face-on profile (C. Peng, 
priv. comm.). However in our case A r so f t > J"b rea k and the de- 
projection is more complex. In order to find a functional form 
for the face-on truncation we built several mock edge-on trun- 
cated disks with GA LFIT and deproj ect the light profiles using 
the "Lucy-method" dWarmelsl 1 19881) . We then fitted these pro- 
files with exponential disks multiplied by the GALFIT trunca- 
tion function. We found that good results are obtain using for 
the face-on view the same A r SO f t as for the edge-on view but a 
flbreak(face-on) ~ '"break + ^A r soft . The effect on the light profile is 
shown in the upper left panel of Fig. [7] 



The bottom panels of Fig. [7] show the rotation curve de- 
compositions for NGC 891 and NGC7814 using the above 
3.6/im-band photometry. We show here the maximum light 
(bulge + disk) solutions for both galaxies. The fit was obtained 
by fitting simultaneously the M/Lssfm ratios of bulge and disks 
to the points inside about 3 scale-lengths (~ 15 kpc for NGC 891 
and ~ 13 k pc for NGC 7814) as it i s customary for maximum 
light tests (Ivan Albada et all Il985l) . The gas surface density 
(middle panels) was multiplied by a factor 1.4 to account for 
Helium. The fitting parameters are given in Table [5] 

In the fit for NGC 891 the bulge dominates in the inner parts 
and accounts for the steep rise and the inner peak of the rotation 
curve. The first point of the rotation curve was included in this fit 
in spite of the uncertainties about its nature and origin pointed 
out above. However, the bulge M/L and the quality of the fit 
would not change significantly, first point included or not. At 
larger radii, beyond ~ 5 kpc, the disk dominates. Its contribution 
to the total mass is about three times as large as that of the bulge. 
The shape of the rotation curve is remarkably well reproduced 
out to the last point. There is no discrepancy between observed 
and model curve and no dark matter halo is required here. This, 
however, is not too surprising. The rotation curve is not very 
extended: it is not traced beyond the bright optical disk and out 
to radii where usually the halo becomes conspicuous. The values 
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Fig. 7. Rotation curve decompositions for NGC 89 1 and NGC 78 14. Top panels: surface brightness profiles at 3 .6 pm for the bulge 
and the stellar disks built from the parameters of Table 2](for NGC 891 the two disks have been added). Middle panels: H I surface 
density. Bottom panels: Best fit without a dark matter halo, the M/L ratios of bulge and disk components are given in Table [5] 
1 kpc = 21.7" and 33.4" for NGC 891 and NGC 78 14 respectively. 



Table 5. Fits to the rotation curves 



Parameter 


NGC 891 


NGC 78 14 


Maximum light 


Bulge M/L 3 . 6(m 
Disk M/L X6/im 
Reduced^ 2 


1.64 + 0.07 
0.90 ± 0.02 
1.05" 


0.64 ± 0.03 
9.25 ± 0.27 
0.59° 


Isothermal halo 


Bulge M/L 3 . 6lml 

Disk M/L 3 , 6llm 

DM halo p (l(r 3 M G per 3 ) 

DM core radius r (kpc) 

Reduced^- 2 


1.63 + 0.25 
0.77 + 0.16 
33.1 ± 16.0 
1.9 + 4.4 
1.30 


0.71 + 0.05 
0.68 ± 0.94 
152.4 + 95.7 
2.1 + 0.6 
0.39 



"Referred to the radial range where the fit is performed (within 3 
scale-lengths), see text. 



obtained for t he M/L ratios o f bulge and disk (Table |5]l seem 
reasonable (cf . I Verheiieni 1 1 997h . 

In conclusion, for NGC 891 a maximum light (bulge+disk) 
solution is quite satisfactory. Bulge and disk together are able to 
explain the rotation curve. Their relative contributions, however, 
are uncertain given the mentioned degeneracy between bulge 
and thick disk. It is clear, at any rate, that in NGC 891 the disk, 
which is responsible for most of the light, contributes also the 
largest part of the mass. 

In the NGC 78 14 maximum-light decomposition the bulge 
dominates in the inner parts. Bulge and disk together provide 
an excellent fit to the rotation curve out to R ~ 13 kpc. Beyond 
that, a mild discrepancy begins to show up between observed 
and model curve. This is the kind of discrepancy found in the 
outer parts of spiral galaxies and usually interpreted as evidence 
for the presence of a dark matter halo. The M/Z^^m for the 
bulge is close to normal, whereas that for the disk is unrealisti- 



cally high. An alternative to such a heavy "dark" disk would be 
a substantial contribution from a dark matter halo (see below). 
The presence of a massive bulge component in this galaxy is 
dictated by the shape of the rotation curve in the inner parts . The 
bulge is completely dominant in the inner 5 kpc. In conclusion, a 
maximum light solution for NGC 7814 provides an excellent fit 
for most of the rotation curve (out to 13 kpc) but it requires an 
unrealistically high M/Z^^m for the disk. A dark matter (DM) 
component is, therefore, needed for the region of the disk and to 
account for the discrepancy in the very outer parts. Clearly, there 
is for NGC 7814 a disk/halo degeneracy. 

3.3. Fits with isothermal DM halos 

The maximum light solutions investigated above are useful to 
understand the possible connection with the light distribution 
and the role played by baryons. They show that the distribu- 
tion of light and that of mass are very similar inside each of 
the two galaxies. Here we add dark matter halos to the fits 
and model the DM distr ibution as a standard isothermal sphere 
(Ivan Albada et all [T985). We fit the four parameters (M/L ra- 
tios and halo parameters) simultaneously. The resulting fits are 
shown in Fig. [8] and the value s of the parameters are g i ven in 
Table Using a NFW profile (iNavarro. Frenk & White! Il997l) 
would not change significantly our results. The best-fit values 
of the concentration parameters are above 10 for both galaxies. 
If we fix them to c = 10 we obtain M/L ratios comparable to 
those in Table [5] except for the M/L of the disk in the free fit of 
NGC 78 14, which becomes M/Z^^m = 2.7; if however we fix 
this latter to 0.7 the fit is still acceptable. 

In NGC 891 the results are very similar to those obtained 
above with the maximum light solution. The fit is still excel- 
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Fig. 8. Rotation curve decompositions for NGC 89 1 and NGC 78 14. Top panels: surface brightness profiles at 3 .6 pm for the bulge 
and the stellar disks built from the parameters of Table 2](for NGC 891 the two disks have been added). Middle panels: H I surface 
density. Bottom panels: Best fit with an isothermal dark matter halo, the parameters of the fits are given in Table[5] 1 kpc = 21.7" 
and 33.4" for NGC 891 and NGC 78 14 respectively. 



lent. The values for the M/Z/^^m of bulge and disk differ only 
slightly and the dark halo plays a minor role. The exclusion of 
the first point of the rotation curve in the fit would not make any 
difference. 

NGC 7814 is more puzzling than NGC 891 because of the 
disk/halo degeneracy. However, a fit with all four parameters 
(M/L ratios and halo parameters) free converges and gives ac- 
ceptable results. The bulge is still dominant in the inner 7-8 kpc 
and very close to maximum with only a minimal change in 
M/Ls.bftm. Further out, the halo now completely dominates in 
the place of the disk. 

These results clearly indicate that the two main luminous 
components in the two galaxies, the disk in NGC 891 and the 
bulge in NGC 7814, are closely linked to the distributions of 
mass as traced by the rotation curves. 

3.4. MOND 

We also compared both rotation curves with th e predictions from 
the Modified Newton Dynamics (MOND) ( iMilgromL Il983h . 
The rotation curve of NGC 891 is reproduced, using the stan- 
dard value of Oq = 3700km 2 s" 2 kpc with best- fit M/L ratios 
for the bulge and disk components of 2.0 and 0.5 respectively. 
The result is shown in the upper panel of Fig. [9] The quality of 
this fit is however not very good leading to a reduced^ 2 of 4.8. 

For NGC 7814 a MOND fit with free M/L 3 . 6flm ratios for 
the bulge and disk components requires for the latter a very 
high value of 4.6. This is needed to fit the outer parts of the 
rotation curve. In this case the M/L3 of the bulge is 0.77. 
If the M/L3 of the disk is kept fixed to a more acceptable 
value o f 1 , the MOND fit using the standard interpolation func- 
tion of lMilgroml (Il983l) is not acceptable (x 1 = 7.2). The bot- 



tom panel of Fig. [9] shows the model predictions for this func- 
tion (thin curve) and also for the so-called "simple" function of 
iFamaev & Binnevl (120051) . Clearly this second function gives a 
much better fit (x z = 1.3 bulge M/L 3 , 6fim = 0.83). We also 
used the "simple" function for NGC 891 but we did not find 
a significant improvement of the fit (disk M/L^.^^m = 0.23, 
bulge M/L 3 6/Jtn = 1.9, x 1 = 4.1). In order to make the fit 
with the standard interpolation function compatible with the data 
of NGC 78 14, the distance to this galaxy should be increased 
by a factor ~ 1.5. This is not compatible with the new error 
determination that allows, a t 3<x, an increase of at most 15% 
dRadburn- Smith et afl 1201 lh . To improve significatly the fit of 
NGC 891 rotation curve the distance should be decreased by a 
factor > 2, incompatible with the new determinations. 



4. Conclusions 

The two galaxies NGC 891 and NGC 7814 are representa- 
tive of two extreme morphologies, with the disk dominat- 
ing in NGC 891 and the bulge almost entirely dominating in 
NGC 78 14. We have derived new rotation curves fo r both . 
Contrary to previo us reports (Ivan der Kruit & Searle 1 Il982l 
Ivan der Kruitl Tl983). the shapes of these curves are found to be 
significantly different. They indicate that in NGC 78 14 the mass 
is more concentrated to the centre as compared to NGC 891. It 
resembles, therefore, the distribution of the luminosity, which 
is more centrally concentrated in NGC 78 14 (bulge) than in 
NGC 891 (disk). 

A decomposition in bulge, disk and halo shows that in 
NGC 891 the disk is the major mass component. The bulge 
contributes about one fourth of the total dynamical mass. In 
NGC 78 14 the bulge dominates almost entirely the total lumi- 
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NGC 891 - MOND 



Fig. 3. Upper: Position-velocity diagram along the major axis 
of the edge-on galaxy NGC 7814. The squares show the rotation 
curves derived separately for the approaching and the receding 
side. The contour levels are: 0.5, 1, 2, 4, 10, 20 mJybeam -1 . 
Lower: Position-velocity diagram along the major axis for the 
model galaxy. 1 arcmin= 4.25 kpc. 
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Fig. 9. Comparison between the rotation curves of NGC 891 and 
NGC 78 14 and the MOND predictions. For NGC 78 14 we show 
the predictions obtained using the standard interpolation func- 
ti on (thin blue curve' ) and th e so-called "simple" (FB) function 
of Famaev & Binnev (2005) (thick blue curve). 



nosity (90 percent of total). In the distribution of mass it is un- 
doubtedly the dominant component in the inner parts. In the 
outer parts the dark matter halo takes over. The disk, unless it 
has an unrealistically high M/L ratio, seems to be a minor com- 
ponent. 

Standard MOND fits do not work perfectly for both galaxies. 
NGC 7814 is well fitted b y the "simple" interpolation function 
dFamaev & BinnevLl2005l) . 

We conclude that in both galaxies, in their bright optical 
parts, the distribution of mass seems to follow closely the dis- 
tribution of light. This implies that either the baryons dominate 
or the dark matter is closely coupled to the luminous component. 
It would be interesting to repeat the same study on galaxies with 
extreme morphologies such as NGC 7814 and NGC 89 1 but seen 
at lower inclination angles. 
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